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ABSTRACT 

Aims. Stellar activity produced by spots and plages affects the radial velocity (RV) signatures. Because even low activity 
stars would produce such a signal, it is crucial to determine how it influences our ability to detect small planetary signals 
such as those produced by Earth- mass planets in the habitable zone (HZ). In a recent paper, we investigated the impact 
of sunlike spots. We aim here to investigate the additional impact of plages. 

Methods. We used the spot and plage properties over a solar cycle to derive the RV that would be observed if the Sun was 
seen edge-on. The RV signal comes from the photometric contribution of spots and plages and from the attenuation of 
the convective blueshift in plages. We analyzed the properties of the RV signal at different activity levels and compared 
it with commonly used activity indicators such as photometry and the Ca index. We also compared it with the signal 
that would be produced by an Earth-mass planet in the HZ. 

Results. We find that the photometric contributions of spots and plages to the RV signal partially balance each other 
out, so that the residual signal is comparable to the spot signal. However, the plage contribution due to the convective 
blueshift attenuation dominates the total signal, with an amplitude over the solar cycle of about 8-10 m/s. Short-term 
variations are also significantly greater than the spot and plage photometric contribution. This contribution is very 
strongly correlated with the Ca index on the long term, which may be a way to distinguish between stellar activity and 
a planet. 

Conclusions. Providing a very good temporal sampling and signal-to-noise ratio, the photometric contribution of plages 
and spots should not prevent detection of Earth-mass planets in the HZ. However, the convection contribution makes 
such a direct detection impossible, unless its effect can be corrected for by methods that still need to be found. We 
show that it is possible to identify the convection contribution if the sensitivity is good enough, for example, by using 
activity indicators. 

Key words. Techniques: radial velocities - Stars: starspots - Stars: planetary systems - Stars: activity 



1. Introduction such an approach compared to the previous ones is that we 

, , , . , 1. 1 , • /, r T-,Tr\ could take the full complexity of the activity pattern into 

^ Among other techniques, the radial velocity (hereafter RV) ^^^^^^^ (structures, spatial distribution, temporal varia- 

, technique has been very successful m detecting exoplanets. ^-^^^^^ ^^^^^^^ ^. ^^^^ ^^^-^^^^ positions 

> , The increasing sensitivity of new instruments (for example ^j^^ ^-^^ partially balance each other out. Our ap- 

■ Espresso on the VLT or Codex on the E-ELT) is expected ^^^^^^ ^^^^ ^^^^^ ^^^^^ temporal behavior of the 

rN to give access to increasingly lower planet masses down • i / u -j-u • j i,- i, i 

' , -2 , X , -,1 1 ,, 1 1 signal (such as with periodograms, which are also very use- 
^ . to Earth-like masses); however, it also makes the observed riJ-iJ-j-Lii u-^r • n j • j 
. . , ' , , , ml tools to detect planets m the RV signal). We derived a 
RV signal sensitive to other sources than planets, such as .^^^ 1 ^^^^^^ ^^^j^ maximum peaks 

stelkr activity and pulsations at ±2 m/s, especially during the solar maximum (R'^x ~- 
The in fluence of spots on th|RV andbisecto^e first ^ ^ . j minimum (R'^^ ^-5.05), the rms 
studied by' Saar fc Donahu^ ^ andEiiii This ^^J^^ ^^/^ ^.^^ me^^tnnm peakf It ±0.6 m/s. Then 
was fol owed by the p recise study o f ti^e influence of a sin- estimated the impact of these spots on the detection of 
gle spot on the signal Desort et al-L M) m various stellar E^^h-like planets and showed that a very good temporal 
conditions such a^jsclmations or angular rotation rates. In ^^^^^ ^^e, ^ iod was mandatory for detecting 
i Lagrange et al. | (M, hereafter Paper I, we simulated the garth-mass planets in the habitable zone, depending on the 
RV signal due to spots that would be detected it the Sun (,ygjg phase 
was observed edge-on over a long period (typically one so- 
lar cycle) with different temporal samplings (between one However, spots are not the only source of RV variations, 
day and 20 days) with only a few gaps. The advantage of The presence of bright plages also induces a variable sig- 

nal in RV, as for spots. They are well-known on the Sun, 

Send offprint requests to: N. Meunier but there are indirect indications that they are also present 
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on other stars, as they sometim es dominate the photome- 
try (e.g. iLockwood et alJ . 120071) . Their contrast in temper- 
ature is less than for spots, but their area is much larger. 
Their contribution is therefore expected to be significant. 
Such a contribution to the RV has never been studied in 
precise detail. Because they are bright, we expect plages 
located at the same position as spots to produce an RV sig- 
nal anti-correlated with the one from the spot (see Paper 
I). However, because they are not localized exactly at the 
same position (in addition to their different spatial exten- 
sion) and have a position-dependent contrast, and because 
in practice there is no one-to-one correspondence between 
the two types of structures (the plage to spot ar e a rati o 
presents a very large dispersion, IChapman et al.l . |2001D . 
the residual signal will be significant. In addition, there 
are many small plage-like structures outside active regions, 
called the network, contributing significantly to the photo- 
metric variations of the Sun. It is crucial to study the two 
contributions (spots and plages) separately and in detail 
because the respective contrib utions of plages and sp ots 
may be different on other stars (jLockwood et al.l . |2007[ ). 

Another contribution of plages comes from the atten- 
uation of the convective blueshift (due to the presence 
of granulation, e.g. iGravl |2009() when magnetic field is 
prese nt. Rough e stimates hav e been mad e on individual 
lines (jSaar fc Don ahue. 1997j, ISaaii l2003l . iKlirster et all . 
l2003f) but the impact on the RV variations derived from 
the cross-correlation function (CCF) has not been studied 
yet. 

The objective of this paper is therefore to estimate the 
contribution of plages to the RV signal (either from the 
photometry or the convection), to compare it to the sig- 
nal produced by spots at the same time, and to study the 
resulting signal. As in Paper I, we consider a Sun seen edge- 
on over a solar cycle. In Sect. 2, we describe in detail the 
data set (spots and plages), the data processing, and the 
observables we obtain. Special care is given to the com- 
parison with the observed photometry (total solar irradi- 
ance) in order to validate the parameters used to make the 
simulation. The resulting RV are analyzed in Sect. 3 and 
compared with other observables (photometry, Ca index) 
in Sect. 4. In Sect. 5, we study the influence of this signal 
on Earth-like planet detection, and conclude in Sect. 6. 

2. Description of the simulations 

In this section we first describe the input data: the list 
of sunspots and plages that is used as input to our sim- 
ulation tool; the way we derive the temperature contrasts 
for these structures; the infiuence of convective inhomo- 
geneities. Then we briefly describe our simulation and the 
produced observables, which are similar to what was done 
in Paper I. 



2.1. Input data 

To get a good temporal sampling of plages, we used 
MDI/SOHO magnetograms (jScherrer et al.L[l995[ ) between 
May 5, 1996 (julian day 2450209) and October 7, 2007 (Ju- 
lian day 2454380). This covers a solar cycle with a tem- 
poral sampling of about 1 day. We therefore used a differ- 
ent sunspot data set from Paper I to cover the full plage 
data set, although there is a significant overlap between 
julian days 2450209 and 2453004. We chose the sunspot 
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Fig. 2. Upper panel: projected spot filling factor versus 
time, in ppm of the solar disk. Lower panel: same for plages. 



grou p data provided by U SAF/NOAA over the same pe- 
riod (http://www.ngdc.noaa.gov/stp/SOLAR/). 

This sunspot data set has a few gaps. Some of them 
correspond to days without any spot while others are true 
gaps in the data. We removed the latter from the MDI data 
set by checking the Wolf number. In the following, there are 
therefore plage structures every day, but in some cases there 
may be no spot on the Sun surface. This leads to 3586 days, 
for a total duration of 4171 days (i.e. a 86% coverage, which 
is similar to what we had for spots in Paper I). 

The spot data are used without any further processing 
and provide 35207 sunspot groups, with the smallest size 
10 ppm (of the solar hemisphere). This size corresponds to a 
radius of 3.1 Mm, assuming a circular shape. We note that 
the surface coverage of the Debrecen data used in Paper 
I is larger by a factor of about 1.57 when comparing the 
same days. The impact of this difference is investigated in 
the next section. 

Plages (in active regions) and network structures (bright 
magnetic structures outside active regions) are extracted 
from MDI magnetograms by applying a threshold of 100 G 
for 9 up to 72° (with 9 the angle between the line-of-sight 
and the normal to the solar surface) in order to avoid the 
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Fig. 1. Upper panel: MDI magnctogram for September 14, 2002 at 4:48 UT [left] and Meudon photospheric spectroheli- 
ogram at 6:41 UT {right). Lower panel: Simulated map showing the plages derived from the magnetogram with a contrast 
increased by a factor 10 for clarity {left) and the spots from the USAF/NOAA data set at 12:40 UT {right). 



noise close to the limb. This threshold is quite conservative 
for limiting the influence of the noi se (about 20 G), an d is 
slightly larger than the one used by iFligge et al.l ()2000D on 
the same data. For the same purpose, the smallest struc- 
tures are eliminated. The smallest size is then 3 ppm (of the 
solar hemisphere), which corresponds to a radius of 1.7 Mm 
assuming a circular shape. This is smaller than the smallest 
spots we consider because we also include network struc- 
tures and not only plages from active regions. This leads to 
1803344 structures, from the quiet network to large plages 
in active regions. The relative contribution of plages and 
network structures is discussed in Sect. 2.3.2. 



The structures extracted from the magnetograms using 
this simple segmentation also include spots of course. With 
the present context, this is not critical, however, since the 
surface coverage of spots is much smaller. Typically, the fill- 
ing factor of spots (determined from the USAF data above) 
represents about 6% of the filling factor of structures deter- 
mined from the magnetograms. Furthermore, modifications 
of the threshold used to identify plages lead to modifica- 
tion of that filling factor by about 10% in the case of large 
structures and up to 20% for small network features, with- 
out changing their visual identification significantly. This 
gives an idea of the uncertainty on the filling factor. The 
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influence of this variation is also investigated in the next 
section. 

The first input to our simulation tool (for both spots and 
plages) then consists in the list of structures: size in ppm 
of the solar hemisphere (before correction for projection ef- 
fects), latitude, and longitude. The projected filling factors 
are shown in Fig. [2] The correlation between the spot and 
plage filling factors is 0.77. This reflects the correlation be- 
tween the two, both on long (cycle) and short time scales 
(rotation), as active regions usually include spots, but also 
the large dispersion of the plage-to-spot size ratio and the 
inclusion of the network structures (not spatially correlated 
with spots as they are located outside active regions). The 
plage filling factor is about one order of magnitude more 
than the spot filling factor. 




1 , , . I , , , I , , , I , , , I , , . : 
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2.2. Contrast choices 

2.2.1. Spot and plage observed contrasts 

There is a large dispersion in the literature concerning 
the spot and plage contrasts. There is also a large tem- 
perature difference from one structure to the next. It 
is therefore difficult to directly use the publishe d con- 
trasts fo und in the literature (for example iFrazieil . Il971l 
iLawrenc e. 1988, ErmoUi ct al., 2007), because they are very 
sensitive to the way plages are defined. The large dis- 
persion between individual plages, as well as the differ- 
ences between plages and n e twork struct u res (iFqukal et aP . 
[T991I iWorden et all [Tool lOrtiz et al.L I2002L I2006D "also 
adds to the uncertainty of these measurements. To real- 
ize a consistent simulation, we computed the photomet- 
ric contribution of each kind of structure (spots, plages, 
network) and compared it with the observed total so- 
lar irradiance (hereafter TSI). In the following we use 
the TSI compile d by G laus Frohlich and Judith Lean 
( |http://www.ngdc.noaa.gov/stp/SOLAR/) fo r the 1996- 
2003 period and described in lFrohlich fc Lead ([TMl)- This 
comparison is therefore made over 2263 points. 

In Paper I, we used a spot temperature deficit of 
550 K, as this lies w ell within the range o f observations 



and s imulations (e.g . Chapman et al.L 11994 iFligge et all . 
I2OOOI IK rivova et al.Ll2003h . The situation is more complex 
for plages. We can either us e a shape such as the one used 
bv lLanza etHI (|2004 I2007D and in Paper I, i.e. a temper- 
ature excess as a function of l-fj, (where /x = cos 9). We 
can also use a slig htly more cor nplex function such as in 
lUnruh et al.) (|1999I) and used by iFhgge et all (|l998l . I2OOOD 
to reconstruct the solar irradiance, to account for the differ- 
ence from zero contrast at disk center. In the following, we 
use a second-degre polyn omial starting with a shape simi- 
lar to the one modeled bv lUnruh et al.l (|l999l) . as it gives a 
slightly better when comparing the simulated irradiance 
to the observed one (see next section). 
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Fig. 3. Upper panel: Plage temperature excess versus fi 
(solid line), and with a modification of the plage size of 
±10% (dashed lines) and ±20% (dotted lines), see text for 
details. Middle panel: Same for the contrast at 600 nm. 
Bottom panel: Plage temperature excess versus the con- 
trast. 



2.2.2. Determination of the spot and plage contrast 

We used the following procedure to determine the spot 
and plage contrasts. We first consider a spot tempera- 
ture deficilEI AT^=-550 K and a plage contrast Cp=0.339- 

^ with respect to the photospheric temperature of 5800 K 
^ defined as {Spi — Sph)/Sph, with Sph the photospheric irra- 
diance and Spi the plage irradiance 



0.563/i ± 0.270/x^ K, w h ich co rresponds to the contrast pro- 
vided bv lUnruh et al.l (|1999[ ) for wavelengths in the range 
470-550 nm. We then computed a contrast at 600 nm as- 
suming a Planck law, which leads to estimate the contribu- 
tion of spots and plages to the TSI. We chose 600 nm for 
reference because that is where the contrast is close to the 
average contrast over 350-1000 n m, representative of the 
bolometric value (|Gondoinl . l2008l ). hence the TSI. 
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Fig. 4. Upper panel: Spot temperature deficit versus tlie 
spot-surface scaling factor for the explored range of scaling 
factors (see text for details). For a given spot-surface scal- 
ing factor, the various points correspond to different plage- 
surface scaling factors. Bottom panel: Same for plages at 
Ai=0. 

This leads to two series, I^p and Ipi, respectively repre- 
senting the contribution of dark and bright features to the 
solar irradiance, in fractions of the TSI of the quiet Sun 
(i.e. contrasts). When multiplied by the quiet Sun irradi- 
ance Sref they allow a direct comparison with the observed 
TSI Sobs (in W/m^) by computing 

Smod — Sref{Isp + Ipl + l)j 

which sums the contributions of spots, plages (including the 
network), and quiet Sun respectively. We can adjust Smod 
to Sobs by applying a scaling factor to each contribution and 
fitting them to minimize 
where 



J2i^mod — Sobs)" / [N — Npar) 



Smod — Sreffrefifsp-^sp ^" fpl-^pl ^" 1); 

and N the number of points, Npar the number of param- 
eters. The parameter fref allows adjustment of the quiet- 
Sun level, while fsp and fpi allow the contrast to be varied 



for spots and plages respectively. For our input spot and 
plage data as described above, we obtain the following best 
choice: S ref= 1365.46 W /m^ (which is very close to the one 
given by iFrohUchl . l2009l for the cycle minimum in 1996), 
ATs=-663 K and ATp=250.9-407.7At-M90.9/.t2 K. Figure [3] 
(solid lines) shows the plage excess temperature as a func- 
tion of /i, as well as the corresponding contrast at 600 nm. 
This contrast is compatible with the literature. The cor- 
relation between the simulated and observed irradiance is 
0.89. 

The final reconstructed irradiance versus the observed 
TSI will be shown in Sect. 2.5, as here we have made a 
simplified computation to correct for projection effects, we 
attribute to each structure the value of jj, corresponding to 
the center of the structure, while in the following we build 
maps of the solar surface. We also neglect the center-to-limb 
darkening (which will be taken considered in Sect. 2.5 where 
we make a more precise computation of Smod)- Figure [S] 
shows the contribution of spots and plages to the irradiance, 
as well as the simulated irradiance used here to derive the 
set of parameters. Our results are quite similar to those 
of lFrohlich fc LeanI (| 19981) . and we reproduce the observed 
irradiance. 



2.2.3. Uncertainties 

To study the influence of our input data set (USAF sunspot 
group and plage extraction from MDI magnetograms), we 
explored the parameter space related to the structure sizes 
by applying a certain scaling factor to them, which should 
represent the uncertainty we have on the structure sizes. For 
spots, we applied scaling factors between 1 (original USAF 
values) and 1.57 (corresponding the Debrecen filling factor), 
i.e. we allow the spot size to vary by about 50%. The smaller 
filling factor in t he USAF catalog may stem from missing 
very small spots (jBalmaceda et al.l . l2009t) . at least partially. 
The correction also included a small dependence on the size, 
as we have been able to compare specific regions, which 
shows that, as expected, the correction is relatively large 
for small structures. For plages, we applied scaling factors 
between 0.8 and 1.2, which takes the assumption concerning 
the inclusion of spots into account, as well as a margin for 
the choice of the threshold. We therefore allowed the plage 
and network sizes to vary by ±20%. 

For each new series of spot and plage sizes, we searched 
for the contrasts that provide the best fit to the observed 
TSI, following the method described in Sect. 2.2.2. Indeed, 
at first order, if one slightly overestimates the plage areas, 
for example, it can be compensated for a larger tempera- 
ture excess. This is true for the reconstitution of the solar 
irradiance and also for RV computations as it also is a pho- 
tometric effect. The quiet-Sun solar irradiance is fixed to 
our best value above, as it is very robust. 

The results are as follows. The corresponding ATg 
ranges between -670 K (for a spot scaling factor of 1) and 
-450 K (for a spot scaling factor of 1.57), with a strong cor- 
relation as expected (see Fig. U). At fi=0, ATp varies be- 
tween ^315 K (for a plage scaling factor of 0.8) and ~215 K 
(for a plage scaling factor of 1.2), and between ~43 K and 
~29 K at /x=l, here again with a strong anticorrelation 
because the areas and temperature excess compensate for 
each other. This anticorrelation is shown in Fig.|4]for /i=0. 
Figure[3]also shows the plage-temperature excess (and con- 
trast at 600 nm) versus fi for the whole domain. Over the 
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whole domain, the remains very low. It does not show 
any systematic effect when changing the plage surface. The 

is, however, smaller for the spot scaling factor of 1 com- 
pared to 1.57. It is also very close to the minimum as a test 
of a few values below 1 for the spot scaling factor shows 
that the then increases again. 

Finally, we note that in Paper I we used the Debrecen 
data set, which corresponds to the spot-surface scaling fac- 
tor of ^1.6 in Fig. |4l We used a spot temperature deficit 
of -550 K, which lies in the middle of our range, instead 
of -450 K. We therefore expect our spot RV signal to be 
weaker by about 20% than that of Paper I, in which the 
spot temperature deficit could not be validated precisely 
using the observed irradiance because plage were not taken 
into account. 

We are therefore confident that our choice of tempera- 
ture contrasts will provide the correct amplitude of the RV 
signal, as we are including two crucial ingredients in our 
simulation: a correct amplitude due to a good combination 
between area and contrast, as checked by the comparison 
with the observed irradiance, and a realistic temporal evo- 
lution of structures thanks to the temporal sampling and 
coverage. 

2.3. Convective inhomogeneities 
2.3.1. Estimation of convective shifts 

In addition to the photometric effect described above, 
another effect can lead to RV and bisector variations. 
The presence of convection in the photosphere produces 
a blueshift of spectral lines, as well as a distort i on of 
these lines (e.e. i Dravins et all . ll98lL iDravinsl . Il982l Il999l 
[Livingston et al. . Il999t ). This comes from the correlation 
between velocities and brightnesses in granules. Granules 
are convective cells at the 1 Mm scale, in which upward 
motions associated to bright area occupy a larger sur- 
face than downward motion area (dark lanes), leading to 
blueshifts and to line distortion. However, where mag- 
netic field is present, it is well know n that convection 
is aff ected, and is greatly reduced (e .g. Hanslmeier et all . 
ll99lLlTitle eTaLl 119871 [Schmidt et al.l . ll988h . This stronglv 
impacts the lines by locally modifying their shape and 
producing a redshift that comp ensates for the convective 
blueshift (ICavallini et al. . 1985L llm merschitt & Schroeter, 
1989llBrandt fc So lankiU 19901. 'Hanslm eier fc Mattia . il99a 
Guenther fc Matt ig. 1991, Hanslmeier et al.l . \1994h . This 



influence of the mag netic fi eld on the lines has also been 
modelled bv .Marguez et al.l (|T996, ) : therefore, we do expect 
this blueshift to be modified by the presence of plages (as 
well as spots, but their contribution is small as they repre- 
sent a much smaller surface and are darker). 

In the following, we estimate the amplitude of the con- 
vective blueshift in the quiet Sun. The convective blueshift 
largely depends on the line depth, and a large dispersion is 
present in the literature as most studies focus on a small 
number of lines. We have used the detailed results of Gray 
(l2009l) for the Sun and derived a blueshift as a function of 



the line depths using a linear function. Then, on the solar 
spectrum obtained by Delbouille, Neven, and Roland (re- 
trieved from the BASS2000 data base), we have identified 
the lines in our spectral range of interest, retrieved their po- 
sition and depth, then simulated a spectrum in which each 
line has the same depth but is shifted by a different amount 
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Fig. 5. First panel: Reconstructed spot contribution to ir- 
radiance. Second panel: Reconstructed plage contribution 
to irradiance. Third panel: Total reconstructed irradiance. 
Fourth panel: Observed irradiance. 



according to its depth and following iGravl (|2009l ). We used 
a linear relation between the line depth and the blueshift, 
which is extrapolated to lines outside the range of depth 
considered byiGray (2009). Finally, we computed the shift 
that would result from a cross-correlation applied on such 
a spectrum. We therefore consider the distribution of line 
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depths in the spectrum and that cross-correlation functions 
(used later to derive the RV) are more sensitive to deep lines 
than weak lines. This produces a convective blueshift cor- 
responding to the integrated light of about 200 m/s that is 
representative of the whole spectrum. The local (vertical) 
convective blueshift is therefore 285 m/s (the correspond- 
ing factor being the result of projection effects combined 
to the c enter-to-limb darken ing). It is similar to th e valu e 
used bv lKiirster et al.l (|2003D derived from Dravins (il999|) , 
but higher than the value considered bv iSaan (j2003i), about 
90 m/s for the vsini we consider. 

iGravl (|2009t ) studied the blue shift only for a qu i et Sun 
(or star). On the other hand, iBrandt fc Solankil (|1990[ ) 
studied the line shape and shifts for different magnetic fill- 
ing factors (including the quiet Sun), by observing at dif- 
ferent positions in a plage. They find a blueshift in plages 
which is 1 /3 of the value in the quiet Sun. We therefore use 
in the following an attenuation of the convec tive blueshitf of 
2/3 of our values, i.e. 190 m/s. The results of lCavallini et al.l 
(fl98l are also compatible with this result. We ex pect the 
attenuation to be around zero in spot umbra (jSolanki 
I1986L iMartinez Pillet et all . I1997D . but the contribution of 
the umbra is very small. Given that the spread in convec- 
tive blueshift associated to plages in the literature, about 
50 m/s, we consider that, in the presence of plages and 
network (see next section), the convective blueshift is at- 
tenuated down to 190 m/s, with an uncertainty of ±50 m/s. 
We neglect in this work the contribution of the Evershed 
effect (mostly associated to spot penumbra) and supergran- 
ulation, because we expect them to only increase slightly 
the jitter on scales from days to weeks. 

2.3.2. Notes on the convective shifts 

All solar observations of the line distortions and shifts due 
to magnetic fields have focused on plages. Our data set in- 
cludes a significant contribution from the network. We also 
know that some abnormal granulation i s associated to it, fo r 
example, to ne twork bright points (e.g . lMuller et al.l . ll989t) . 
More recentlv. iMorinaea et al.l (|2008 ) have shown the sup- 
pression of convection around small magnetic concentra- 
tions; however, there is no specific study of the lineshift 
associated to the network that we are aware of. In the fol- 
lowing, we consider the RV signal associated to all struc- 
tures, but will also estimate it for regions of a signifiant size 
only in order to test the possible range of variation. There 
is of course a continuum in size from the largest activ e re- 
gions to the small network structures (jMeunie 3. 1200311 . but 
only considering sizes above 100 ppm (corresponding to a 
radius of 10 Mm) would probably give a good estimation 
of a lower limit for the convection contribution to the RV. 
This is discussed in Sect. 3.3.2. 

Finally, it should be noted that the final RV signal due 
to convective inhomogeneities is sensitive to the actual ar- 
eas of structures and to the AV we consider. As before, 
for area and temperatures, there is a trade-off between the 
two. However, unlike the temperature, we do not have the 
possibility of testing it independently. 

2.4. Simulated spectra 

Simulated spectra are computed as described in Paper I. 
A rigid rotation of 1.9 km/s at the equator is assumed. 




1364 1365 1366 1367 1368 
Observed irradiance (W/m^) 

Fig. 6. Total reconstructed irradiance versus the observed 
irradiance (W/m^). 



For each type of structure (spots and plages), we build a 
3D hemisphere map of the visible Sun seen edge-on. We 
assumed circular shapes both cases, mostly for simplicity 
as we do not expect the shape to significantly impact on 
our result, and more importantly to be able to extrapolate 
easily to other stellar cases in future works. For the pho- 
tometric contribution of magnetic structures, the resulting 
spectrum is then computed assuming that spots (respec- 
tively plages) emit a solar-like spectrum following a black 
body law with a temperature as described in the previous 
section. To take the convection into account, we consider a 
solar-like spectrum (black body with the photospheric tem- 
perature) shifted by a velocity perpendicular to the solar 
surface (which produces a redshift) of 190 m/s. 

2.5. Computation of the observables 

The o bservables are then computed using the SAFIR soft- 
ware r^all and et all . [2005) as on actual stellar spectra. In 
this paper we focus on the analysis of the RV. As in Paper 
I, computations were made only on one spectral order (or- 
der #31) as the whole process is very time consuming. 
We therefore obtained three series of RV, due to spots, to 
plages (and network) through their photometric contribu- 
tion (hereafter the plage signal), and to the infiuence of 
plages on convection (hereafter the convection signal) . The 
bisector velocity span (hereafter BVS) due to convective 
inhomogeneities has a very complex influence on the final 
line shape, and is not characterized with sufficient precision 
in the literature. We therefore only analyze the RV in this 
paper. 

We also computed the photometry at 600 nm to check 
that we obtain the expected agreement. Figure El shows the 
good correlation (0.89) between the final simulated pho- 
tometry and the observed TSI. The residual difference be- 
tween the model and the simulation largely comes from 
the strong temperature variation from one structure to the 
next. However, we capture the essential variability of the 
photometry. 
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Table 1. RV rms (in m/s) for specific components. 



Earth-like planets detection capabilities. 





spots 


plages 


sp+pl 


conv. 


total 


All 


0.34 


0.31 


0.33 


2.38 


2.40 


Low 


0.09 


0.10 


0.08 


0.44 


0.44 


High 


0.48 


0.44 


0.42 


1.39 


1.42 



Table 2. RV rms and peak-to-peak amplitude (in m/s) and 
relative rms photometry for the three periods. 
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2.6. Simulation of an active star with a planet 

To test the planet detectability providing the activity we 
model, we added the signal of a 1 Mearth planet located at 
1.2 AU (see Paper I) to the activity signal, as well as some 
random noise at various level. This is analyzed in Sect. 5. 

3. Results 

In this section we present the different contributions to the 
RV signals over the whole period. In addition, in order to 
have a better idea of the dependence on the activity level, 
we also show the results for two 9-month periods as in Paper 
I: a period of low activity from July 1, 1996 (JD 2450266) 
to April 1, 1997 (JD 2450540), leading to 213 data points 
over 274 days, and a period of high activity from February 
1, 2000 (JD 2451576) to November 1, 2000 (JD 2451850), 
leading to 274 data points, also over 274 days. They will 
be denoted as "Low" and "High" , respectively, in the next 
figures and tables. 

In addition to the total RV and the 3 individual com- 
ponents, we also discuss the contributions due to spots 
and plages separately (the total being denominated as the 
spot-f plage signal) and compare them with the convection 
signal. This is important as they have a different origin 
(even if the convection signal relates to plages): on a star 
with similar plages but much weaker convection, the plage 
signal would remain the same but the convection signal 
would become very small. 

3.1. RV variations 

3.1.1. The spot and plage induced RV variations 

We first consider the individual time series of RV due to 
the photometric contributions of spots and plages. They 
are shown in Fig. [T] The amplitude of the spot signal is 
similar to what we found in Paper I (slightly smaller as ex- 
pected), with maximum RV up to a few m/s during high 
activity periods. The RV is lower than 0.6 m/s during the 
minimum. The plage signal is, as expected, also higher dur- 
ing the high activity period (up to 1.6 m/s), i.e. about a 
factor two below the spot signal. The plage signal during 
the low activity period reaches 0.4 m/s, slightly below that 
of spots. It should be noted, however, that during that pe- 
riod the spot signal may be zero for an extended period of 
time, while the plage signal is always different from zero. 
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Fig. 8. Upper panel: RV signal from spots (dotted line), 
plages (dashed line), and spots-|-plages (solid line) in m/s, 
for Julian days 2451790-2451850, i.e. during the high ac- 
tivity period. Lower panel: Same for spots-f plages (dotted 
line), convection (dashed line), and total (solid line). 



Figure |8] (upper panel) shows the different contributions in 
more detail over a short period (60 days) . As expected, the 
spot signal is partially compensated for that of plages (it 
was already true in the very simple computation made in 
Paper I), but not entirely. This can be explained as most 
active regions contain spots (leading to a compensation) 
but the area ra 1 :io be tween spots and plages varies (e.g. 
IChapman et al.l . |2001[ ) and there are plages with no spot. 
Over the whole data set, the correlation between the spot 
and plage RV is indeed -0.47. 

Figure ini shows the distributions of RV for the spot and 
plage signals. The distributions arc not Gaussian, as they 
show extended tails, and, in the case of spots, a sharp peak 
at zero RV (related to the significant number of days with 
no spot). The shape of the distribution is related to the 
shape of the solar cycle, i.e. to the variation in the activity 
level over time. A similar size distribution of structures with 
a constant activity level gives a Gaussian distribution. A 
sinusoidal shape for the solar cycle, for example, only leads 
to a small departure from the Gaussian shape while the 
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Fig. 9. Histogram of RV (m/s) for spots+plages (solid line), 
spots (dotted line), and plages (dashed line), for the whole 
period. 



proper shape gives a distribution similar to what is shown in 
Fig.m The rms RV for the whole period and the two periods 
corresponding to low and high activity are summarized in 
Table [TJ The sum of the spot and plage RV is indicated 
as "sp+pl". It confirms that at all time spots and plages 
provide similar rms, but as they partially compensate for 
each other, the resulting rms is in fact similar to that of the 
individual components. 

Figure [TU] shows the periodogram for the RV signal over 
the whole period for the spot+plage signal. The peaks cor- 
responding to half the rotation period are strongly empha- 
sized compared to the peaks at the rotation period. It may 
stem from the presence of active longitudes separated by 
180° as discussed in Paper I. 

3.1.2. The convection-induced RV variations 

The convection signal is much stronger than the spot and 
plage ones as shown in Fig. [7] and Fig. HI It is also funda- 
mentally different, as it is always positive so that all con- 
tributions add up. Therefore plages at different longitudes 
do not compensante for each other (as in the previous case) 
and the signal is maximum when the features are at the cen- 
tral meridian. This signal is then correlated very well with 
the filling factor of plages (correlation close to 1). When 
considering a small time series, however, as in Fig. [51 even 
if the signal is shifted because of the convection signal, the 
rms is not very different from the one caused by spots and 
plages alone (only a factor 2 larger this time) . The influence 
on the observed RV will therefore depend on the temporal 
coverage compared to the length of the cycle, as shown in 
Fig. [TT] As previously. Table [T] shows the rms for the RV 
due to convection. While considering restricted periods of 9 
months (low and high activity periods), the rms of the con- 
vection RV remains significantly larger than the one caused 
by spots and plages. We note a factor 6-7 between the rms 
due to spots and plages and the rms due to convection when 
computed over the whole series or at cycle minimum. The 
factor is twice smaller during cycle maximum. 
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Fig. 10. Upper panel: Periodograms for the spot-|-plage RV 
for the whole serie. The horizontal dotted line show the level 
for a false- alarm probability of 1%. Lower panel: Same for 
the total RV, i.e. taking the inhibition of the convective 
blueshift into account. 



3.1.3. The total RV variations 

Finally, we consider the total signal as shown in Fig. [71 
The corresponding rms and maximum amplitude (peak-to- 
peak) are shown in Table [H for the whole period, as well 
as the low and high activity periods. For the whole period, 
we reach an rms of 2.4 m/s, with a peak-to-peak amplitude 
above 11 m/s, clearly dominated by the convection signal. 
The amplitude is significantly lower during the low activity 
period, but can still reach peak-to-peak amplitudes of 2 m/s 
with a rms around 0.4-0.5 m/s. 

Figure ITU] shows the periodogram for the total RV signal 
over the whole period. The convection RV, which dominates 
the total signal, emphasizes the rotation period, however, 
because the sign is always the same. For the same reason, 
the long time-scale features are also emphasized. As the 
rotation period is often not known precisely, any more than 
the zero on the RV curves, such a difference between the 
two kinds of contributions may be difficult to separate. 
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Fig. 11. Histogram of RV (m/s) for convection, for the 
whole period (solid line), low activity period (dotted line), 
and high activity period (dashed line). 

3.2. Comparison with observed solar RV 

There are no observations of the solar RV during the period 
covering our data set (cycle 23) that we are aware of. We 
therefore considered three previous studies covering either 
cycle 21 or 22. We have already compared the RV signal 
caused by spots in Paper I. It was not possible to use these 
observations to validate our simulation. By including plages 
it may be possible to make a more precise comparison, how- 
ever^ 

iMcMillan et all ()1993D measured the solar RV using the 
moon light observed in the violet part of the spectrum. 
They have determined an upper limit for intrinsic solar vari- 
ations of ±4 m/s, during cycle 22 essentially (1987-1992). 
This is less than our total RV variation. However, it is quite 
compatible with our rms RV for the spots and plages only 
(0.33 m/s), as well as our maximum peak-to-peak ampli- 
tude (2.4 m/s), even taking the slightly larger contrasts 
in the violet into account. In the wavelength domain they 
used, ther e is a very large proportion of deep lines, but as 
shown by iGravl (2009) , these are the lines that have the 
lower convective blueshift, so it is possible that their small 
RV comes from their being predominantly sensitive to the 
spot-f plage RV. 

During a similar period (1983-1992), 

iDeming fc: Plvmatd ()1994[ ) measured the solar RV in 
the 2.3 /Ltm domain. They find a long-term variation with 
a peak-to-peak amplitude of about 30 m/s, and a large 
dispersion on short time scales (10-20 m/s). Considering 
that cycle 22 was more active than cycle 23 by about 30%, 
we would expect long-term variations of about 11 m/s, i.e. 
3 times less than their observation. It is quite possible that, 
in the wavelength domain they consider, the convective 
blueshift is much more pronounced, but there has been no 
such study that we are aware of to confirm that possibility. 
Their short-term variations are also significantly greater 
than our results for spots and plages, especially considering 
the reduced contrast expected at 2.3 /im. 

Finally, Jimenez et al. ( 1986|) measured the solar RV in 
the potassium line at 769.9 nm (|Brookes et al.l . I1978D dur- 
ing cycle 21 (1976-1984). On long time scales, they find 



a peak-to-pe ak amplitude of 30 ni/s, w hich is compatible 
with that of iDeming fc Plvmatd (|l994f ) and also 3 times 
larger than our convective signal. Their short time-scale 
variations can reach 20 m/s during cycle maximum, and 
are therefore also significantly more than our results. 

3.3. Discussion on the convection RV amplitude 

3.3.1. Impact of the convective blueshift amplitude 

We computed the convective RV variations on a short pe- 
riod (60 days) when using a AV from 50 m/s to 300 m/s. 
It appears that both the RV and BVS amplitudes vary lin- 
early with AV. This means that it is easy to extrapolate 
to different convection conditions from our results. 

3.3.2. Impact of the areas cancelling the convective blueshift 

As pointed out in Sect. 2.3.2, the reduction of the convec- 
tive blueshift in the network structures outside active re- 
gions is more uncertain than in plages. We therefore recom- 
pute the convection RV while selecting structures (plages) 
only above 100 ppm. This corresponds to structures above 
305 Mm^. We recall that the convective RV is correlated 
very well to the filling factor, and therefore we find results 
that are consistent with the relative cont ribution of thes e 
structures over the solar cycle. As shown bv lMeunie 3 (l2no3h . 
all structures including the smallest one (at the MDI spa- 
tial resolution) are correlated with the solar cycle, and the 
amplitude of variation decreases when the size decreases. 
Here, we find that these structures contribute to 35% of 
the rms RV during the low activity period, and to 55 % 
during the high activity period. 

For comparison, the same large structures contribute in 
a larger proportion to the plage (photometric) RV: 68 % 
during the low activity period and 83 % during the high 
activity period. The reason is that, in that case, because 
similar structures on both side of the central meridian can- 
cel each other and because small structures are more nu- 
merous and more homogeneously distributed over the solar 
surface, small structures contribute less than large struc- 
tures, relative to their contribution to the filling factor. For 
convection however, all contributions add up, so they are 
directly related to the filling factor. 

We also find that very small structures, for example 
those smalle r than 80 Mm^ (typical network structures for 
IWand . ll988l) . should contribute to 30 % of the filhng factor 
at cycle maximum and 70 % at cycle minimum. When con- 
sidering structures only above that size, we find that the 
plage RV has almost the same rms as before (0.29 m/s in- 
stead of 0.31 m/s, confirming the insignificant contribution 
of the smallest strructures to the plage signal) , but the rms 
of the convective RV is about 2 m/s instead of 2.5 m/s, so 
that these same small structures significantly contribute to 
that signal. As pointed out above, this only estimates the 
various contributions, but the exact attenuation of the con- 
vective blueshift in these regions is not well-known, even if 
we know that the granulation associated to them is abnor- 
mal. 

Finally, the conservative threshold we have used to 
extract plage and network structures from MDI magne- 
tograms i mplies that we are missing a significant part of the 
fiux (e.g. iKrivova fc Solanki [200I in very quiet regions, 
which may also be associated to a small attenuation of the 
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Fig. 12. Total RV (m/s) versus the observed TSI (W/m^) 
for the whole data set (dots), for the low activity period 
(yellow triangles) and for the high activity period (pink 
squares) . 



Fig. 13. Total RV (m/s) versus the Ca index for the whole 
data set (dots), for the low activity period (yellow triangles) 
and for the high activity period (pink squares). 



convective blueshift, with a badly defined level. Therefore 
it is likely that our RV signal slightly underestimates the 
exact value. 

In conclusion, considering plages above 100 ppm only 
provides a conservative lower limit for the convection con- 
tribution, i.e. about half the total signal during cycle max- 
imum and about a third during the minimum. 

3.3.3. Impact of the method 

To separate the photometric contribution and the RV due 
to the reduction of the convective blueshift in plages, we 
computed these two contributions separately. We also com- 
puted the RV signal on a short time series for plages that 
take both effects into account simultaneously; i.e., we con- 
sider a plage contribution shifted by Ay and with the usual 
plage contrast in temperature. We want to compare the re- 
sulting RV: RVaii with the sum of the plage and convection 
RV computed before RVsum- The difference between the two 
present an rms of 0.04 m/s, which is one order of magnitude 
smaller than the plage RV for that period (about 0.4 m/s) 
and of course much smaller than the convection contribu- 
tion. RVaii and RVsum are very well correlated (correlation 
of 0.97). There is however a systematic offset of 0.23 m/s 
between the two, with RVaii larger than RVsum- This value 
is, however, stiU smaller than the average RVaii (5.8 m/s). 
We conclude that our approach leads to the correct ampli- 
tude and variability of both contributions. 

4. Search for correlations of the radial velocity 
signal with other observables 

4.1. Total solar irradiance 

We compared the obtained RV with the observed TSI over 
the period 1996-2003. As photometry is often used as an 
activity criteria, it is important to check how the RV re- 
lated to activity correlates with the TSI. There is no cor- 
relation between the RV due to spots and plages and the 



TSI. Low TSI values tend to correspond to periods with a 
small rms in RV, as already pointed out by Kiirster et al.l 
()2003D . When considering the total RV (Fig.[ni), there is a 
small correlation between the RV and TSI over the whole 
period (correlation of 0.40), although this correlation tends 
to disappear when considering shorter periods. 

4.2. Ca index 

Another commonly used indicator of stellar activity is the 
Ca index. It is therefore useful to study the correlation be- 
tween this indicator and the observed RV. As in Paper 
I we use the Ca index provided by the Sacramento Peak 
Observatory, for the period 1996-2002. Their is no corre- 
lation between the RV signal considering only spots and 
plages and the Ca index. However, as already noted in 
Paper I for spots, periods with a low Ca index correspond 
to a small RV dispersion, while periods with larger Ca index 
correspond to a larger RV dispersion. When considering the 
total RV, i.e. taking the convection into account, over the 
whole period (Fig.[T31), the correlation is strong (0.89). This 
is not surprising as the Ca index is related to plages, and 
the total RV is dominated by the convection signal (also 
related to plages, with no sign change). 



5. Planets 

To test the detectability of earth-type planets in the hab- 
itable zone, we have considered, as in Paper I, a 1 Moarth 
planet orbiting at a 1.2 AU from the star on a circular or- 
bit. The RV amplitude of such a planet, for a system seen 
edge-on, is 0.08 m/s and its orbital period 480 days. We 
computed the resulting RV signal for the whole data set, 
as well as for two periods of low and high activity. The 
noise level i s chosen to be 1 cm/s (the goal for Espresso 
on the VLT, iD'Odorico fc the CODEX/ESPRESSO te^ . 
2003) ] 5 cm/s, and 10 cm/s (for Codex on the E- 
ELT, iP'Odorico fc the CODEX/ESPRESSO teaml . [2001 . 
Different temporal samplings were considered. 
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Fig. 14. Upper panel: Periodogram of the spot+plage RV 
(m/s) added to the planet (see text) and a 1 cm/s noise, 
over the whole period. The sampling is the original one 
with 4-month gaps every year. The vertical dashed line in- 
dicates the planet orbital period. The dotted lines indicate 
the power corresponding to false alarm probabilities of 1% 
(upper line) and 10% (lower line). Lower panel: Same for 
the total RV added to the planet (see text) and a 1 cm/s 
noise. 



5.1. Spot and plage RV variations 

Figure [T3] first shows the periodogram of the signal corre- 
sponding only to spots and plages, with a 1 Mcaith planet 
orbiting at a 1.2 AU and a 1 cm/s noise level. The peak cor- 
responding to the planet is highly significant (false alarm 
probability below 10~^^). The sampling is the original one 
with 4-month gaps every year to simulate a star that would 
not be visible at all times (dec~ —45°) as in Paper I. The 
same is true when considering only a few years of low activ- 
ity (extended to cover 790 days) or high activity (extended 
to cover 2000 days). When degrading the temporal sam- 
pling down to 4 or 8 days, the amplitude of the planet peak 
decreases and becomes barely significant, especially on the 
shorter periods, as in Paper I. 



5.2. Total RV variations 

When considering the total RV (including the convection 
contribution), again with the 1 Mcarth planet orbiting at a 
1.2 AU, a 1 cm/s noise level and a 4-month gap every year. 
Fig. [HI shows that the peak corresponding to the planet is 
not significant. It is also adjacent to larger peaks at nearby 
periods. This is due to the fact that convection dominates 
the spot-|-plage contribution, the noise, and the planet, so 
that the total RV shows a long-term variability related to 
the solar cycle, including on the scale of the year or a few 
years, hence the adjacent peaks (see for example Fig. [TU|) . 

If we decrease the amplitude of the convection contri- 
bution by just a factor two (see Sect. 3.3.2), the planet 
peak has a false alarm probabiblity of 18%; for a factor 5 
it becomes 3 10^^. The latter is very significant, but it is 
still adjacent to nearby peaks that are as significant, there- 
fore it would be very difficult to identify it as caused by a 
planet. Increasing the noise level or degrading the tempo- 
ral sampling does not change the nature of the problem. In 
conclusion, we would need to decrease the convection con- 
tribution by one order of magnitude to be able to identify 
the planet signal, even when using the best sampling. 



6. Conclusion 

We computed the expected solar RV variations caused by 
the photometric contribution of spots, as well as the con- 
tributions of plages through their photometry and the sup- 
pression of the convective blueshift between 1996 and 2007. 
The Sun was considered to be seen edge-on and observed 
in conditions similar to HARPS (spectral coverage and res- 
olution). This work complements Paper I, which only con- 
sidered the spot contribution. Our simulation reproduces 
the observed photometry, so we are confident that their 
relative contributions are correct. Our approach takes the 
complex activity pattern observed on the Sun into account. 
Even if the full characterization of such a complex pattern 
is currently beyond reach for other stars, precise photo- 
metric observations by COROT shows that solar-like stars 
can exhi bit several short-liv ed spots at the same time (for 
example iMosser et al.l . l2009l) or show a complex pattern be- 
tween spots and plages fe.g. lLanza et al.l . l2009l for a young 
solar- like star). We obtain the main following results. 

The photometric contributions of spots and plages only 
partially compensate for each other, so that the sum of 
these two signals is not very different from the spot signal. 
The shape of the RV distribution is not Gaussian and may 
provide information on the shape of the cycle, which may 
be a useful diagnosis tool in the future. 

The plage contribution due to the convective blueshift 
suppression dominates the total signal. Unlike the previous 
one, the sign is always the same (a redshift when activ- 
ity is present), i.e. is additive, and its long-term amplitude 
is about 8 m/s. The short time-scale RV rms is close to 
the m/s (and down to 0.4 m/s during the low activity pe- 
riod). Finally, we point out that this con tribution is very 
line- dependent, as shown bv iGravl (|2009[ ). and therefore a 
precise computation must be made depending on the spec- 
tral range considered. If many lines are available, it should 
be possible in principle to compare the RV series computed 
using the strong and weak lines separately in order to iden- 
tify the origin of the variation. 
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The comparison with observations of the solar RV is 
not easy, because of different temporal coverages and wave- 
length domains. The photometric signal should be sensitive 
to the wavelength domain because the contrast (assuming 
a black body) is lower in the red than in the violet part of 
the spectrum. On the other hand, the convective blueshift 
is very line-sensitive. Here we have estimated the contribu- 
tion when considering most lines of the visible spectrum. 
On the other hand, all solar RV observations have focused 
on a single line or a restricted part of the spectrum, for 
which this convective blueshift might be very different and, 
unfortunately, has been poorly studied. 

The RV is not correlated well with the total solar irradi- 
ance. Therefore a lack of correlation between observed RV 
and photometry should not lead to the conclusion that the 
observed RV cannot be due to activity. This is true for both 
the photometric contribution of spots and plages to the RV 
and for the convection signal. On the other hand, when 
the total RV is dominated by the convection signal, we ex- 
pect from the present results a strong correlation with the 
Ca index. It is also dependent on the variability of the Ca 
index. The correlation remains very small, however, when 
only considering the photometric contribution of spots and 
plages to the RV. 

We also tested the influence of activity and noise on the 
detection of an earth-like planet at 1.2 AU on a circular 
orbit. When considering only the photometric contribution 
of the spots and plages, the results are similar to that of 
Paper I. A good temporal sampling is necessary to detect 
the planet. On the other hand, when also taking the con- 
vection into account, the activity- related RVs completely 
dominate the signal at a level of about 8 m/s, and it seems 
impossible to detect the planet. Because there are other sig- 
nificant peaks in the same range of periods, the convection 
signal would need to be at least one order of magnitude 
weaker to make the identification of the planet peak possi- 
ble, otherwise the adjacent peaks would still be of similar 
amplitude. 

In conclusion, the photometric contribution to RV of 
both spots and plages at a level comparable to the Sun 
should not prevent detection of earth-like planets in the 
habitable zone around solar-like stars, providing a very 
good temporal sampling (better than 8 days and over a long 
monitoring period) and a good signal-to-noise ratio. On the 
other hand, the attenuation of the convective blueshift due 
to active regions, at the same level, seems to make it im- 
possible to detect such a planet in similar conditions, as 
the orbital period is in the same range as the long-term 
strong activity signal. However, the properties of this con- 
tribution are such that some diagnosis may be established 
more easily than with the photometric contributions. 

— The convective blueshift presents a strong spectral sen- 
sitivity, so that using two masks (of respectively deep 
and weak lines) should lead to significantly different RV 
variations. For example, for a AV of 200 m/s over all 
lines, we find 455 m/s for lines with depths between 0.5 
and 0.95 (weak lines) and 155 m/s for lines with depths 
below 0.5 (strong lines), i.e. an expected factor 2.7 be- 
tween two such RV series. 

— The correlation with the Ca index is very strong. Such a 
correlation could be used to correct long -term RV vari- 
ations, as done bv ISaar fc Fischeil (|2000f ). 



— The BVS remains to be studied in detail, but we can 
expect it to be of large amplitude given the impact of 
that contribution on the lines, allowing it to be above 
the noise of future instruments. This will have to be con- 
firmed using a realistic modelization of the line distor- 
tion when taking the convective blueshift into account. 

All these diagnosis will be very sensitive to the noise level. 
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Fig. 7. Left column: RV (m/s) for spot signal, plage signal, convection signal, and total (From top to bottom), for 
whole data set. Middle column: Same for the low activity period. Right column: Same for the high activity period. 



